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Abstract 

With the availability of large-scale redshift survey 
data, it is now becoming possible to explore correla- 
tions between large-scale structure and the properties 
and morphologies of galaxy clusters. We investigate 
the spatial distributions of a 98% complete, volume- 
limited sample of nearby (z < 0.1) Abell clusters with 
well-determined redshifts and find that cooling flow 
clusters with high mass-accretion rates have nearest 
neighbors which are much closer than those of other 
clusters in the sample (at the 99.8% confidence level), 
and reside in more crowded environments out to 60 
h^Q Mpc. Several possible explanations of this effect 
are discussed. 

Subject headings: cooling flows — galaxies: clusters: gen- 
eral — large-scale structure of universe 

1. Introduction 

Intriguing correlations between cluster properties 
and larger-scale structures have been noted over the 
past two decades. For example, the major axes of 
Abell clusters may be aligned with nearby clusters 
(e.g., Bingelli 1982; Plionis 1994), X-ray substruc- 
ture within clusters appears to point to neighboring 
Abell clusters (West, Jones, & Forman 1995), and 
the tails of Wide- Angle Tailed (WAT) radio sources 
within clusters tend to be aligned with the long axis of 
the nearest supercluster (Novikov et al. 1999). All of 
these alignment effects are broadly consistent with the 
scenario of hierarchical structure formation whereby 
clusters form through the flow of material along sheets 
and filaments (e.g., Shandarin & Klypin 1984; Pauls 
& Melott 1995; Bond, Kofman, & Pogosyan 1996; 
Colberg et al. 1999). Here we investigate another 
possible link between cluster properties and larger- 
scale structures: the environments within which cool- 
ing flow clusters are located. 

Cluster cooling flows are believed to result when 
gas within the central ~ 100 kpc cools radiatively 



on a time-scale less than a Hubble time; the resul- 
tant decrease in pressure allows gas to flow towards 
the core (see, e.g. the review of Fabian 1994). Sev- 
eral groups have estimated that (60 — 90)% of clusters 
contain cooling flows (Edge, Stewart, & Fabian 1992; 
White, Jones & Forman 1997; Peres et al. 1998) with 
mass deposition rates (obtained by detailed deprojec- 
tion analysis of their X-ray images) ranging from a 
few to more than one thousand M Q yr _1 . The en- 
vironmental factors which determine whether or not 
a given cluster will develop a cooling flow are still 
unclear and have not been properly addressed by nu- 
merical simulations. However, the development of a 
strong cooling flow was observed in the single cluster 
simulation of Katz & White (1993), and detailed ID 
simulations incorporating many realistic effects sug- 
gest that strong cooling flows are a natural outcome of 
cluster formation over a wide range in masses (Knight 
& Ponman 1997). 

Though cooling flows may be a generic feature of 
cluster formation, it appears that the recent merger 
history of a cluster is closely linked to the existence 
and/or strength of a cooling flow. McGlynn & Fabian 
(1984) and Fabian & Daines (1991) suggested that 
cooling flows may be disrupted by cluster mergers. 
This is supported by the fact that there are no cool- 
ing flows in several massive clusters which appear to 
have undergone recent mergers (e.g., Coma: Burns 
et al. 1994; A2255: Burns et al. 1995). Unfortu- 
nately, the observational picture is not entirely clear 
as there are also cooling flow clusters which appear 
to have undergone a merger (e.g., A1664: Allen et 
al. 1996; A2597: Sarazin et al. 1995). Nonetheless, 
the observations are consistent with a scenario in 
which recent and/or massive mergers tend to disrupt 
any pre-existing cooling flow. For example, Buote & 
Tsai (1996) demonstrated that cooling flow strength 
is anti-correlated with the degree of X-ray substruc- 
ture (which is likely a measure of recent merger activ- 
ity). Recent numerical simulations (Burns et al. 1997; 
Gomez et al. 1999) of cooling flow mergers indicate 

that high-Af flows are more difficult to disrupt than 
low-M flows and that the likelihood of disrupting a 
given cooling flow increases with the mass of the in- 
falling cluster. 

To date, there has been no systematic investigation 
of the large-scale environments within which cooling 
flow clusters are found. A possible connection be- 
tween large-scale structure and cluster cooling flows 
comes from the evidence that mergers can disrupt 
cooling flows which, in turn, suggests that massive 
cooling flows have been relatively undisturbed for 
some time. As Abell clusters are thought to delineate 
the supercluster distribution rather well (e.g. Bah- 
call et al. 1994), we decided to investigate the spatial 
distribution of Abell clusters around cooling flows. 

2. Procedure 
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2.1. The Abell Cluster Sample (AC) 



We began by constructing a complete sample of 
all rich Abell/ACO clusters (Abell 1958; Abell, Cor- 
win, & Olowin 1989) with R > 1, m 10 < 17.0, 
|6| > 30°, and spectroscopically-determined redshifts 
0.012 < z < 0.10 (where R is Abell's richness class 
and mio is the magnitude of the tenth brightest 
galaxy within the cluster). The cut in galactic lati- 
tude rejects regions of the sky where obscuration and 
poor sampling may plague the ACO catalog while our 
magnitude limit corresponds to z <~ 0.13 assuming 
an mio-z relationship (e.g., Batuski & Burns 1985). 
The total number of ACO clusters satisfying these 
cuts is 284, of which we kept the 277 with measured 
redshifts. The majority of the cluster redshifts come 
from the ESO Nearby Abell Cluster Survey (Katgert 
et al. 1996) and the MX Northern Abell Cluster Sur- 
vey (Slinglcnd et al. 1998), although <~ 50 are unpub- 
lished redshifts from Miller et al. (1999b). 

It is important to note the level of completeness 
(98%) and quality of this Abell/ACO dataset. Miller 
et al. (1999a) and Peacock & West (1992) have shown 
that R > 1 Abell/ACO cluster suffer much less 
from spurious cluster selection (due to line-of-sight 
anisotropics) than do samples containing R = clus- 
ters. By limiting the samples to R > 1 clusters we 
are sacrificing some poorer clusters with published 
mass deposition rates. However, since we are inter- 
ested in nearest neighbors, it is imperative that our 
base sample be significantly complete in redshift mea- 
surements. A magnitude-limited and volume-limited 
subset of R > 1 Abell/ACO clusters is the only vi- 
able dataset for an analysis such as this. This sam- 
ple of Abell/ACO clusters has a nearly constant spa- 
tial number density (n = 9 x 10~ 7 /)| Mpc~ 3 ) out to 
z-0.10. 

After selecting an initial sample based on the above 
criteria, we removed any clusters which were closer 
to a boundary of our volume than to their nearest 
neighbor in order to carry out the nearest-neighbor 
analysis described below. The remaining set of 202 
clusters, which shall be referred to as AC, is a 98% 
complete, volume-limited sample of Abell clusters. 

2.2. The Cooling Flow Clusters 

Using NED, we searched the astronomical litera- 
ture for papers dealing with AC clusters in order to 
find published determinations of M. The single most 
comprehensive source of clusters with detailed X-ray 
deprojections and analysis is the 207 clusters with 
Einstein data studied by White, Jones, & Forman 
(1997; hereafter WJF). Our second major source is 
the analysis of ROSAT data performed by Peres et 
al. (1998) on a sample of the 55 X-ray brightest clus- 
ters originally identified by Edge et al. (1990). For 
clusters that appear in both the WJF and Peres et 
al. lists, we use the ROSAT-derived values of M from 
Peres et al. (taking the PSPC-derived value when- 
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Table 1: The CF (cooling flow) sample of clusters. 
These are all the AC clusters with non-zero pub- 
lished values of M. The references for the quoted 
values are: P (Peres et al. 1998), S (Sarazin 1986), 
and W (WJF97). The quoted values assume H o =50 
km/s/Mpc. 
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ever available). A handful of additional M values 
were gleaned from other papers (Sarazin 1986; Edge 
& Stewart 1991; Pierre & Starck 1998). 

All told, 55 of our AC clusters have published val- 
ues of M and 26 of these have non-zero mass ac- 
cretion rates (see Table [[]). The fact that half our 

clusters with estimated M have non-zero mass ac- 
cretion rates is consistent with the cooling flow oc- 
curence rate estimated by WJF. For the analysis dis- 
cussed below, we created several different subsamples 
of cooling-flow clusters; CF1 contains only those clus- 
ters with M > 50M Q yr _1 , CF2 is a sample with 
M > 35M yr _1 , and CF3 is a sample of those clus- 
ters with M — (jj^j > 20MQyr _1 (where av, is the 
estimated l-a lower uncertainty). 

2.3. The Non-Cooling Flow Clusters 

The WJF catalog quotes M — for 27 of our 
AC clusters. However, after taking into account the 
spatial resolution of the images, WJF conclude that 
only 10 of these are truly excluded as being cooling 
flows. We take these to be our most secure sample of 
non-cooling flow clusters (NCF1). This sample com- 
prises: A0019, A1185, A1213, A1291, A1656, A1913, 
A2040, A2147, A3158 and A3744. Our second sample 
of non-cooling flow clusters, NCF2, contains all those 

AC clusters in WJF which have M = and an up- 
per uncertainty less than 20. After removing A2065 

which has M / according to Peres et al. (1998), the 
NCF2 sample is made up of: A0150, A0154, A0168, 
A0399, A0500, A0690, A1185, A1213, A1377, A1656, 
A1809, A1913, A2040, A2079, A2092, A2124, A2147, 
and A3744. 

2.4. Selection Effects 

Our AC sample is well-suited for the spatial distri- 
bution analyses we are interested in. Unfortunately, 
the 55 AC clusters with published values of M form 
an inhomogeneous sample since they were observed 
and studied by various authors for different reasons. 
Nonetheless, we find (§3) that subsamples of cooling 
flow and non-cooling flow clusters drawn from this 

sample of clusters with known- M differ significantly 
in their nearest-neighbor distributions. It is difficult 
to imagine how such an effect could arise spuriously if 
the clusters with known-M form a random sampling 
of the AC clusters. We point out that the clusters 
with known-M cannot be selected from AC on the ba- 
sis of an X-ray flux or luminosity cut suggesting that 
they are more-or-less randomly selected. Though all 
the AC clusters are potential members of the X-ray 
Brightest Abell Clusters (XBACs) catalog (Ebeling et 
al. 1996), only 61 AC clusters and only 37/55 of the 
AC clusters with known-M meet the XBACs flux- 
limit (even our two most massive flows do not make 



the catalog). 

3. Results 

We wished to determine if the existence of cool- 
ing flows is influenced by the proximity of other clus- 
ters. We tabulated the distance to the nearest neigh- 
bor for every cluster in our AC sample and then 
considered whether the distribution of these nearest- 
neighbor distances is statistically different for any of 
our subsamples. Since clusters are correlated, we ex- 
pect these distances to be smaller than the 104/i^q 1 
Mpc mean distance between our clusters. Distances 
to all clusters were determined assuming go = and 
H Q = 50 km s" 1 Mpc" 1 . 

Since we are looking for a systematic difference 
in the typical distance of the nearest neighbor, the 
Wilcoxon Rank-Sum test is appropriate (Lehmann & 
D'Abrera 1998). This test consists of putting all the 
nearest- neighbor distances in rank order. The sum of 
the ranks of each subset is examined for consistency 
with the null hypothesis that there is no difference 
in the rankings. We emphasize that the significance 
given by this test explicitly accounts for fluctuations 
that may arise due to sample size. 

Each subsample was compared with its own con- 
trol sample made up of all the AC clusters not in 
that particular subsample. The results of the analysis 
are given in Table |2| including the significance of the 
result that clusters in the sample are systematically 
closer to their nearest neighbor than those in their 
control group (AC minus the subsample). Although 
it is not used in the assessment of significance, we also 
tabulate the mean distance to the nearest-neighbor 
(nnd) for each subsample and its corresponding con- 
trol sample. The analysis shows one very significant 
result; massive cooling flows (CF1) tend to have much 
nearer closest neighbors than their control clusters. 
In fact, the average nearest-neighbor distance (nnd) 
is just half what it is for the remainder of the AC 
sample. When less massive cooling flows are included 
in the subsample (CF2 and CF3), the effect becomes 
less pronounced (though still significant at the 95% 
level). The significance is lowest for CF2 which, by 

definition, includes some clusters whose M values are 
consistent with zero. We performed the same analysis 
with the non-cooling flow samples to make sure that 
this is not a spurious effect. Indeed, the most secure 
sample of non-cooling flow clusters (NCF1) is not sig- 
nificantly different from the remainder of the AC sam- 
ple. The second non-cooling flow sample (NCF2) con- 
tains some potential cooling flows (since M = Otjf ) 
and has a higher likelihood of differing from its control 
sample than NCF1. This suggests a possible trend in 
which nnd decreases with increasing M. 

As a double-check, we calculated the mean num- 
ber density of clusters in spheres of specified radii 
centered on the clusters of the entire AC sample and 
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Table 2: Results of the Wilcoxon Rank-Sum test com- 
paring nearest-neighbor distances (nnd) of our vari- 
ous subsamples with the rest of the AC sample. S 
is the significance (%) that the subsample is system- 
atically more crowded than the control. The average 
nearest-neighbor distance (nnd) for both the subsam- 
ples and the control samples are also tabulated. 



on those in the CFl sample. For any given radius r, 
only those AC (or CFl) clusters which were at least 
a distance r from a boundary of the AC sample vol- 
ume were considered. We used all 277 ACO clusters 
to find density inside a given radius. The results are 
depicted in the top panel of Figure [j] and show that lo- 
cal mean densities (for r < 601b io 1 Mpc) are higher for 
the CFl clusters (dashed line) than for the entire AC 
sample (solid line) . The significance of the difference 
is assessed in the bottom panel which shows the prob- 
ability (using a Kolmogorov-Smirnov test; Lchmann 
& DAbrera 1998) that the distribution of densities 
at a specific radius could result from a random se- 
lection of the AC parent population. The differences 
are highly significant out to 60 h^ 1 Mpc, although 
roughly what would be contributed by an average of 
one extra neighbor near each CFl cluster. 

Our main result is the close nearest neighbor and 
generally increased density of the large-scale envi- 
ronment for large distances around massive cooling 
flow clusters. One might worry that this effect could 
be a richness bias as our CFl cluster sample has a 
larger fraction of R > 1 clusters than the AC sam- 
ple, causing the CFl sample to be more strongly 
spatially correlated (Kaiser 1984) simply on this ba- 
sis. After removing the CFl clusters, the AC sample 
contains 155 R=l clusters, 33 R=2 clusters and 2 
with R=3 whereas the CFl sample has 6 R=l and 
6 R=2 clusters. We therefore constructed four "pro- 
portioned" subsamples of AC that each contained all 
33 of the R=2 clusters as well as 33 randomly cho- 
sen R=l clusters. Averaging our Wilcoxon Rank- 
Sum test for these four samples produced a 99.5% 
confidence that the CFl sample is more crowded 
than these richness-proportioned control groups. The 
mean nearest neighbor distance for the proportioned 
controls is (34-38) Mpc as opposed to 18.4/iJq 1 
Mpc for CFl. This, along with the diagnostics for 
the proportioned set in Figure 1 indicate that the CFl 
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Fig. 1. — The top panel shows the mean number 
density of galaxy clusters in spheres of specified radii 
centered on the clusters in: the AC sample (solid 
line), the CFl sample (dashed line), and the propor- 
tioned AC sample (dotted). More significantly, the 
bottom panel shows the probability that the densi- 
ties around CFl clusters (whose mean is shown in the 
upper panel) could result from a random selection of 
the AC sample (solid line) or the proportioned AC 
sample (dashed line). 



sample is significantly different from both AC and a 
richness-proportioned subsample of AC. We conclude 
that any richness bias in our result is a small effect. 

4. Conclusions 

We find that the most massive cooling flows (those 
with M > 50M Q yr _1 ) have significantly closer nearest- 
neighbors (nnd = 18.4/i^q 1 Mpc) than the typical 
cluster (nnd=38 Ah^Q Mpc) in our complete, volume- 
limited sample of rich Abell clusters. We cannot ex- 
clude that this result is related to some unknown se- 
lection effect connected with why the clusters were 
originally selected for X-ray observation. However, 
this seems unlikely given that Table 2 shows a trend 
for the significance of the result to decrease as less 
massive cooling flows arc added to the subsample and 
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to disappear completely when samples of non-cooling 
flows are used. 

How can one interpret our result that massive cool- 
ing flow clusters tend to reside in environments that 
are crowded with other clusters? One might ex- 
pect that mergers would be much more common in 
such regions and that therefore cooling flows would 
be less likely to survive (e.g. Fabian 1994). How- 
ever, it is possible that these regions are crowded be- 
cause they have not yet finished collapsing and their 
member clusters have not recently experienced vio- 
lent mergers. By this interpretation, many clusters 
which appear isolated may have recently gobbled up 
their neighbors and should show signs of dynamical 
activity. 

Since cooling times will decrease as density in- 
creases (at fixed temperature), we might expect to 
find cooling flows in regions of high baryon density 
which would typically be regions of high matter and 
cluster density (Bardeen et al. 1986). However, cool- 
ing flows are found not only in clusters, but also in 
anemic small groups and individual galaxies. 

Cooling flow clusters are usually highly symmetric 
indicating dynamical relaxation. Thus one might ex- 
pect cooling flow clusters to be old clusters, formed 
from high-amplitude perturbations which collapsed 
early. Especially with the kind of cosmological power 
spectra that appear viable at this time, one would 
expect a high-amplitude peak to be surrounded by a 
larger region of high density. 

We note that the density of clusters around a cool- 
ing flow cluster is significantly higher out to about 60 
h^gMpc. This is interestingly close to the lengthscale 
for supercluster effects on WATs as found in Novikov 
et al. (1999) and to the wavelength scale going non- 
linear today. It is suggestive of a general link between 
cooling flows and the current collapse of larger-scale 
perturbations. It may be related to a coherent gas 
flow into the clusters on an even larger scale than 
previously suspected. 
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